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Abstract. Motivated by recent discovery of strongly spin-orbit coupled two- 
dimensional (2D) electron gas near the surface of Rashba semiconductors BiTeX 
(X=C1, Br, I), we calculate thermoelectric responses of spin polarization in 2D Rashba 
model. By self-consistently determining the energy-and subband-dependent transport 
time, we present an exact solution of the linearized Boltzmann equation for elastic 
scattering. Using this solution, we find a non-Edelstein electric-field induced spin 
polarization which is linear in the Fermi energy Ep, when Ep lies below the band 
crossing point. The spin polarization efficiency, which is the electric-field induced spin 
polarization divided by the driven electric current, increases for smaller Ep. It is 
shown that, as a function oi Ep., the temperature-gradient induced spin polarization 
continuously increases to a saturation value when Ep downs below the band crossing 
point. As the temperature tends to zero, the temperature-gradient induced spin 
polarization vanishes. 

Keywords: thermoelectric response, Rashba spin-orbit coupling, Boltzmann equation, 
analytical solution 


1. Introduction 

Two dimensional (2D) electron systems with spin-orbit coupling (SOC) show a great deal 
of fascinating transport phenomena due to the mixing of the spin and orbital degrees of 
freedom, providing the possibility of realizing all-electrical and all-thermal spin control 
in semiconductor structures. These are the main topics of the rapid developing research 
helds of spintronics [T] and spin-caloritronics j2]. In the electrical spin control, the 
generation of a spin current and a nonequilibrium spin polarization transverse to an 
applied electric held without external magnetic held are remarkable, known as the spin 
Hall ehect [3lll] and electric-held induced spin polarization [5H7], respectively. Their 
thermal counterparts toward all-thermal spin control, i.e., the spin Nernst ehect [8HI2] 
and temperature-gradient induced spin polarization HMS], have also attracted more 
and more interests recently. 

The 2D electron system (2DES) with Rashba SOC has been one of the most widely 
used models to investigate aforementioned ehects [El UHl IISI [H] . In the 2D Rashba 
model, two bands cross at zero energy, one of them is always positive and the other one 
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possesses a band valley regime below the band crossing point as shown in Fig. 1. In this 
valley regime the dispersion curve is not monotonic in momentum space. This regime 
possesses nontrivial topology of the constant energy surfaces (or Fermi surfaces) [TH] . 
which leads to some exciting theoretical predictions, e.g., the enhanced superconducting 
critical temperature [I6], the non-Dyakonov-Perel spin relaxation behavior im and the 
signihcantly enhanced room-temperature thermoelectric hgure of merit [18] . There have 
been a few theoretical studies [m[T71 - l20] on the transport properties when the Fermi 
energy is in or near the band valley regime. However, in Rashba systems formed in 
conventional narrow-gap semiconductor heterostructures [21], the Rashba spin splitting 
energy is so small that the band valley structure can not survive the weak disorder 
broadening and thermal smearing even at very low temperatures. In these systems the 
Fermi energy usually lies quite above the band crossing point, therefore the band valley 
is irrelevant to transport. 

Recently experimental progress has been made by the discovery of giant bulk 
and surface Rashba SOC effects in V-VI-VII polar semiconductors BiTeX, (X=C1, Br, 
I) [ 22 H 23 ]. In these noncentrosymmetric semiconductors, hrst-principles calculations 
and ARPES measurements have clearly demonstrated the existence of 2DES conhned 
near the surface with giant Rashba energy as large as about lO^meV [25l[26]. In such 
2DES, the investigation of electrical and thermal spin control is of signihcance due to the 
giant SOC which is promising for spintronics and spin-caloritronics applications. While, 
when the electron-impurity scattering dominates, for the case that Fermi energies lie 
below or in the vicinity of the band crossing point, the relaxation time approximation 
(RTA) used in previous theoretical works on the nonequilibrium spin polarization [mi20] 
may not work well due to the giant Rashba SOC. This motivates us to systematically 
investigate the thermoelectric response of spin polarization in 2DES with giant Rashba 
SOC, focusing on the consequences of different Fermi surface topologies between the 
two sides of the band crossing point. 

In this paper, we employ the semiclassical Boltzmann equation to calculate the 
spin polarization induced by electric held and temperature gradient. We focus on the 
2D Rashba model at low temperatures where the static impurity scattering dominates. 
Our calculation is based on an exact transport time solution of the Boltzmann equation 
in the Born approximation, different from the widely used modihed RTA and constant 
RTA schemes [27]. We show that the electric-held induced spin polarization (EISP) as 
a function of the Fermi energy Ep behaves diherently between the two sides of the band 
crossing point Ep = 0. A linear dependence of EISP on Ep is obtained for Ep < 0, 
dihering from the Edelstein result [5] for Ep > 0. The spin polarization efficiency, 
dehned as the ratio between the EISP and the driven electric-current, increases for 
lower Ep. The temperature-gradient induced spin polarization (TISP) is calculated, 
and its dependence on the Fermi energy, changing from large positive values to be quite 
below the band crossing point, is continuous and monotonic. It is also shown that the 
temperature-gradient induced spin polarization tends to zero at vanishing temperatures. 
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2. Semiclassical Boltzmann descriptions of thermoelectric spin responses in 
Rashba 2DES 


2.1. Basic solutions for the 2D Rashba model 


We study the 2D Rashba model with spin independent disorder 

^ ■ (P X z) + (r), (1) 

where R (r) = Ei(r - R*) is the disorder potential produced by randomly 
distributed ^-scatters at R* and is assumed to be standard white-noise disorder: 
(Ek'kl^)^jg = Here Him is the impurity concentration, 14;'k the spin-independent 

part of the disorder matrix element and (-Eis the disorder average, m is the in-plane 
effective mass of the conduction electron, p = hk the momentum, a = cry,az) are 
the Pauli matrices, a the Rashba coefficient. Eigenenergies of the pure system are 
Exk = + Aak, with inner eigenstates |MAk) = [1, —exp where A = ± 

and tan0 = ky/k^. The wave number at a given energy E > 0 in the A band is given 
as kx{E) = -Xkn + ^xM. + 2EjiE (see Fig. 1), where we define the Rashba energy 
Er = m and kn = The density of state (DOS) at a given E > 0 is given by 
N>{E) = ZxNx{E) where 


Nx (E) = No- 


kx (E) 


’kx{E) + Xkn 

Here Nq = is the DOS of 2D spin polarized parabolic spectrum. For F > 0, the 
group velocity and intraband spin matrix element are given by 

No hkx{E) 


v(F, A,0) = 


Nx{E) 


m 


(3) 


and 


{uxkx{E)W\uxk^iE )) = A (sin0x- cos0y), (4) 

respectively. 

The direction of the group velocity is the same as the corresponding momentum 
due to the isotropic band structure and the monotonic E — k curve when E > 0. The 
directions of spin in the two bands at the same polar angle 0 are opposite to each other. 

The lower band has a valley centered at k^, and the DOS has a one-dimensional 
(ID) character in the {kn) < E < 0 regime [16] with (kn) = —\Er the 
energy of the bottom of dispersion curve. For (fc^j) < F < 0 there are two wave 
numbers k -2 {E) < kn < fc_i {E) with k^i, (E) = kji + (— ^x/Ef + 2EEji, where 
—u = — 1, —2 denotes the two monotonic branches in this energy regime (see Fig. 1). 
The DOS A^< {E) in the band valley regime is given by A^< {E) = Ev=i N-i, {E) where 

k-u {E) 


N-AE) = N,y 

\k-u [E) — ku] 

For E_ {kn) < F < 0, one obtains 

{u-k-XE)WW-k.^{E)) = - sin 0X-F COS 0y 


(5) 


( 6 ) 
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Figure 1. Band structure of the 2D Rashba system. The energy of the band crossing 
point is chosen to be zero. The wave number and energy of the bottom of the dispersion 
curve is and respectively. Corresponding to a given energy E > 0, the wave 

number in ± band is denoted by k± (E). For —I^Er < E < 0, there are two monotonic 
regimes on E — k curve: the one from k = 0 to k^ is marked by the branch —2, 
whereas the other from k = k^ to 2k n marked by branch —1. The wave number 
k-i, (E) represents the wave number in the —v branch at given E, where ^' = 1, 2. 


and 


v(E,= (-1)^ ^ 


No hk, (E) 


(7) 


N_i, (E) m 


The direction of the gronp velocity is parallel (anti-parallel) to the corresponding 
momentnm in — z/ = — 1 (—2) branch, respectively. This is cansed by the non-monotonic 
E — k curve in the band valley regime. The directions of spin in the two monotonic 
branches at the same (j) are the same. These characters show different spin-and group 


velocity-textures of constant-energy circles between the band valley regime and the 


> 0 regime. 

Exactly at the band crossing point {E,k) = (0,0), the eigenstate as well as the 
group velocity and spin matrix element are not well-dehned because the polar angle 0 
is arbitrary. However, this does not bring any influence on physical quantities since the 
DOS at this point is zero, as shown in Eqs. ([2]) and (|5]). As for {E,k) = (0,2/cjj), the 
group velocity and spin matrix element are both well-dehned and continuous. 

2.2. Basic formulas of nonequilibrium spin polarization and Boltzmann equation 

The out-of-equilibrium spin density response to external helds can be obtained in 
the semiclassical version of linear response theory, which can be decomposed into the 
intrinsic and extrinsic parts 


(<^) = : 


( 8 ) 


where 


d)ini = 



( 9 ) 
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Here is the equilibrium Fermi-Dirac distribution function (DF) and gi denotes the 
out-of-equilibrium change of DF, I = (A, k) is the eigenstate index denoting the band 
and momentum. |'^®) is the eigenstate of the disorder-free Hamiltonian in the absence 
of external fields, \Sifi) describes the virtual interband transition induced by the weak 
external helds. 

originates from the intrinsic mechanism based solely on the spin-orbit conpled 
band structure. It is not difficult to verify that = 0 for the present model, so one 

only needs to analyze the extrinsic spin density response which depends on the 

existence of disorder via gi. 

gi can be calculated by the semiclassical Boltzmann eqnation in the presence of an 
nniform weak electric held and small gradients of chemical potential and temperature 
in nonequilibrium steady states. Here we consider low temperatures where the static 
impurity scattering dominates the electron relaxation. The Boltzmann eqnation reads 
dfO 

=-'^wu^iigi-gi,], ( 10 ) 

^ ii 

where the generalized force acting on the state / is F; = — V/r -|- eE with E, 

/i, T being the electric held, chemical potential and absolute temperature, respectively. 
wi,j is the transition rate from state /' to /, which can be determined by the golden rnle 
in the qnantum mechanical scattering theory. In the present system withont anomalous 
Hall ehect, the lowest order Born approximation is snfficient: 

wu,i = \{ui,\ui)f 6 {El - El,) , (11) 

ToNq 

with To = ^ No ^ ^ When E > 0, the intraband and interband elastic scattering 

can be represented by coy’f {E = Ei) = J dEiiWup. 

(E) = [1 + AA' cos (0' - 0)] . (12) 

Whereas for E_ {kjf) < E < 0 we introduce {E = Ei) = J dEiiWy^i to represent 

the intra-branch and inter-branch scattering: 

= ;^^[l + cos(0'-0)]. (13) 

3. The exact solution of the Boltzmann equation 

In this section we analytically solve the Boltzmann eqnation based on the isotropic 
transport times. For E > 0, the Boltzmann eqnation includes both direct intraband 
and interband elastic scattering; while for E_ {kn) < E < 0, only intraband scattering 
in the lower band occurs. Due to the band valley structnre below the band crossing 
point, the solution in this regime is nontrivial and completely different from ordinary 
single-band cases. Finally we clearly show that for positive and negative energies, the 
DFs are formally similar. 
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3.1. The exact solution of the Boltzmann equation for E > D 
When > 0, the Boltzmann equation can be re-expressed as 

FE-y(E,\,^)aEf = -Y."yiE) I(E) 

X [g, {E, tl (v (E, A, 0))) - gx> {E, ^ (v (E, A', 0')))] , (14) 

where -3 (v) denotes the angle of the direction of v with respect to that of the applied 
generalized force = —-^^VT—Vp-l-eE. Due to Eq. ([3]), "d (v {E, A, 0)) = -3 (Ica {E)). 
Above Boltzmann equation can be solved by introducing the isotropic transport time 
for electrons with energy E in the A band as 

^A {E, (kA (E))) = {-dEf) F^; ■ V (E, A, 0) Ta (E) , (15) 

and the transport time is determined self-consistently by substituting Eq. flTSl) into Eq. 
(ITill . Thus we obtain 




E (E) 


A' 


2ti 


(E) 


1 — cos (0' — 0) 


TV (E) 

d(e) 


(16) 


where we use the relation |v (E, A', 0')| = |v (E, A, 0)| suitable for Rashba 2DES. In fact, 
Eq. flT6|) contains two coupled linear equations determining r+ and r_, yields 

2 


Ta (E) = To 


Nx (E) / 2iV, 


N. 


iV>(E) 


where iV> (E) = 2Nq. Combining with Eq. 
compact form 

jj(B,^(kj(B))) = (-SE/“) Fe 


( 17 ) 

[3]), the out-of-equilibrium DF takes a 
hkx (E) 


m 


-T-o, 


which satishes the particle number conservation requirement 

d(j) 


E j dEN, (E) J 


2ir 


gA(B,!l(kA(B))) =0. 


(18) 


(19) 


Eq. flTSl) looks similar to the DF for spin degenerate free electron gas without SOC 
but the important difference is that in the present case the group velocity is given by Eq. 
(]3|) rather than , for positive energies. This solution is the same as that obtained by 
employing a custom-designed ansatz for the DF [29] in anisotropic Rashba-Dresselhaus 
2DES, while our approach is based on the simple physical picture of isotropic transport 
time on constant-energy circles. On the other hand, this transport time solution is 
different from the modihed relaxation time approximation (MRTA) solution for isotropic 
multiband systems: where 


1 


MRTA 


T, 


= 


1 - 


|Vf/| 

|v/| 


cos {-d (v0 - -d (v;,)) 


( 20 ) 


Comparing Eq. (120]) with (ITOil . it is obvious that this MRTA solution can not be 
self-consistently obtained from the Boltzmann equation for Rashba 2DEG, because if 
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one substitutes above into the Boltzmann equation, it is Eq. fllbp rather than 

Eq. fl20l) that will be arrived at for Thus the MRTA only makes sense as an 

approximate solution for Rashba 2DES. We can further point out that, in the present 
model, the result of Eq. fl20|) is the same as Eq. (fT7|) only in the zeroth order of SOC 
while different from the latter even in the 1st order of SOC in the small SOC limit. 
Thus this MRTA in Rashba 2DEG even can not be regarded as a better one than the 
constant RTA obtained by neglecting the scattering-in term directly. 


3.2. The exact solution of the Boltzmann equation for {ku) < E < 0 

When E_ [kfi) < E < 0, hj converting the momentum integration in Eq. (fTOj) into 
energy integration and noticing the different orientations of group velocity in the two 
monotonic branches, the Boltzmann equation can be re-expressed as 

Fe ■ V {E, S SeP = 

v' 

X [g_^ (E, -d (v (E, -z/, 0))) - g_^, (E, -d (v (E, -z/', 0')))] , (21) 

which is similar to Eq. flT^ for E > 0. The derivation of the transport time solution of 
Eq. fl^ is thus similar to that when E > 0. Substituting 

g., {E, (v (E, -iz, 0))) = {-dEf) Fe ■ v (E, -iz, 0) r_, (E) (22) 

into Eq. fl^T]) . taking into account the fact that in the band valley regime the direction 
of the group velocity can be parallel or anti-parallel to that of the momentum, i.e., 
■d (v {E, —1, 0)) = 13 (k_i (E)), i3 (v (E, —2, 0)) = i3 (k _2 (E)) + vr and then 


cos'd (v {E,-iz',(j)')) _ cos^ (k_^.' (E)) 

cos'd (v (E,—z/, 0)) cos'd (k_i. (E)) ’ 

we get the following self-consistent equation for 


( 23 ) 


r-. (E) 


N-.' (E) j 


d(j)' 


P-T, [E) 


X 


1 -(-!)"■" cos (0'-0) 


T-u' {E) 


r-. {E) 

Here we have used the relation |v (E, — z/', 0')| = |v (E, — z/, ( 
is found as 

2 


(E) = To 


N_,{E) f 2W 


No \N^iE) 


(24) 

Then the transport time 

(25) 


where ( 


N<{E) 


_ Ef+2EiiE 


. Comparing this transport time for negative energies to 
Eq. (ITTD for positive energies, one can see that they share the same form. 

Therefore the nonequilibrium DF satisfying the particle number conservation 
requirement is 


(E) = (-Se/“) F 


E • 


-1) 


hk^ {E) 


m 




2No 


N<{E) 


(26) 
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Here and below we use the simplified notation gx (E) and {E) to represent the DF 
for brevity. It is obvious that this DF for negative energies is formally similar to that 
for positive enegies when the latter is re-expressed as 

2No 


gx (E) — {—dEf ) Fe •- To 


m 


N>{E) 


(27) 


except one significant difference: the (—1)^”^ factor for negative energies. This factor 
denotes nothing but the important fact that, for electrons with negative energy on the 
branch z/ = 2, the group velocity and momentum has the opposite directions (see Eq. 
([7]) and Fig.l). 

By Eqs. (1T5]) and (12B]) . the out-of-equilibrium DFs in above two energy regimes are 
continuous at E = 0: ( 7 + {E —)■ 0+) = g_2 {E —)■ 0“) = 0, g_i {E —)■ 0“) = g_ {E ^ 0+). 


4. Electric-field and temperature-gradient induced spin polarization 

Since the nonequilibrium state is driven by the effective electric field E* = E — 
and temperature gradient (—VT), the spin density response takes the following form in 
the linear response regime 

{S')ext = Xe-E*+ XS7T ■ . (28) 

Here the EISP coefficient Xe and TISP coefficient Xs 7 t can be calculated from the 
2nd equation of Eq. ([9]), where the momentum integration is performed by integrating 
over energy and polar angle. Substituting the out of equilibrium DF, i.e., Eqs. (1T8|1 
and (126|) . the Xe and Xvr containing contributions from both bands are given by 
Xe {T, fi) = Xe,+ {T, /i) + Xe,- {T, /i): 

XE.+ = dEN+ (B) (-9e/“) V (B, +, ^) 

(E) ('Uk+{£:)|o-|Mk+(E)), 

Xe.- (U fi) = e y ^ r (-S) i-SEf) V (B, 

xr_ (E) {uk_(E)W\uk_{E)) 

+e[^f dEY,N-AE)(-aEn ( 29 ) 

J 27r JE_^kn) , 

XV (E, -I/, (j)) (E) (nk_4E)|cr|Mk_.(E)), 

and 

1 poo rp _ 

XVT (T. A*) = - / dE {-dsf) -Ef^XE [E). (30) 

Here and below we use the simplified notation Xe {E) to represent the zero-temperature 
EISP coefficient Xe (T = 0, E) for brevity. Xe is a tensor and has two indices: Xe (hi) 
where i specifies the spin component, and j the direction of the electric field. Due to the 
isotropy, we can apply the generalized force only in x direction for the calculation, and 
only XEiipx) will be calculated below (xe(^,^) = 0 ). Also we will drop the indices 
{y,x) in Xe (y,£), for simplicity, xe- 
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EfIEr 

Figure 2. Spin polarization efficiency. 


4 . 1 . EISP 

The zero-temperature EISP for Ep > 0 can be obtained easily from Eq. fl29|) as 

CY 

Xe {Ep > 0) = (31) 

This result has been well-known since Edelstein |5]. It is independent on the Fermi 
energy, because the directions of spin on the inner (-I-) and outer (—) Fermi circles are 
opposite at the same polar angle 0 and the E^r-dependence of EISPs in both Fermi 
circles cancels. 

While, for Fermi energies below the band crossing point, the EISP takes the 
following non-Edelstein form 

Xe {Ep < 0) = eTo—2No ■ (2^) 

It is linearly dependent on the Fermi energy, different from the positive Fermi energy 
case. In the band valley, the orientations of spin on the inner (—2) and outer (—1) 
Fermi circles are parallel at the same 0 and the Fir-dependence of EISPs of both Fermi 
circles does not cancel. Since the Fermi surface topology in the band valley differs from 
that above the band crossing point, the behaviors of EISPs are different between the 
two regimes, xe (F) is continuous at F = 0. The contribution to xe (0) entirely comes 
from the outer Fermi circle {Ep = 0,k = 2kp), since the DOS at the band crossing point 
{Ep = 0, A; = 0) is zero. 

We compare Eqs. fl3lll and fl32l) to the EISP obtained by employing the constant 
RTA [H]: {Ef > 0) = erf 2iVo, Xe^^ {Ef < 0) = erpNo^lE^, with r the 

constant relaxation time independent on the energy and band. For F^r > 0 the EISP 
obtained by the constant RTA has the same Edelstein form, while when Ep < 0 the 
constant RTA result shows different FiT’-dependence from Eq. fl32ll . 

Now we calculate the spin polarization efficiency, dehned as the ratio between 
the electrical-held induced spin polarization and the driven electric-current density. 
The electrical conductivity for the same model has been given by [2H]: u{Ep > 0) = 

+ Therefore, the spin 

polarization efficiency is given by a single expression suitable for both positive and 
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negative Fermi energies 
Xe {Ep) 


hk 


R 


<t{Ef) eBsl + li:’ 


(33) 


which increases for decreased Fermi energies (shown in Fig. 2). Eqs. (El]), fl3^ and 
show that when the Fermi energy lies below the band crossing point, alough the EISP 
is lowered, higher spin polarization efficiency is achieved. 


4.2. TISP 

We substitute xe (E) into Eq. 




ksT 


= X, 


P 


ksT 


= —t 


, and define 
Ep 


knT 


— —to 


a{ti) = 


dx 




dx 


x,b{ti) = 


dx 


'h 


dx 


X 


then the TISP is found as 


_ kp , . 2kBT TT^ 

XvT - —Xe (0) — - — 

e Er 3 


1 - 


b{ti) -tia(ti) 
TT^/S 


(34) 


(35) 


The relation between t 2 and ti (i.e., the chemical potential at low temperatures) 
can be obtained by the consideration about electron density 


t 2 

t 2 


ti — 0, Ep ^ kpT, 
'ksT 


U = o 


t2 —tl — 


71 


Er 
kpT 


, \Ep\ ~ fcsT, 


(36) 


6 Er + 2Ef 


—Ep 3> kpT. 


Here we only consider small thermal fluctuations kpT <C Er and Ep + \Er S> kpT, so 
that the band valley structure and the Fermi surfaces survive the thermal smearing. In 
some new materials with giant Rashba effect, e.g., strongly spin-orbit coupled 2DES near 
the surface of Rashba semiconductors BiTeX (X=C1, Br, I), Er is about 35 ~ 200meR. 
Therefore, kpT Er and Er + ^^Er 3> kpT can be satisfied at low temperatures about 
several Kelvins for not too low Fermi energies. 

Consequently we can set ^2 = H in the expression for xwt since Eq. fl35ll has already 
beeno(^): 


a TT 

XVT = To-2NokB — 
n 3 


1-3 


b (^ 2 ) - t2a (^ 2 ) 


71^ 


2 kBT 
E 


R 


(37) 


This is our main analytical result for TISP. This formula is still valid for Fermi energies 
near the band crossing point, where the Sommerfeld expansion is not suitable for the 
treatment of Eq. fEOl) . due to the fact that the energy-derivative of Xe{E) is not 
continuous at the band crossing point. 

According to Eq. fE7|) . xvt is given in Fig. 3, in units of Tof2NokB^^^^-. When 
Ep/kpT > 5, xvT almost vanishes. When Ep/kpT < -5, xvT/rof2NokB^^^ 
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Ef/ZcsT 


Figure 3. Temperature-gradient induced spin polarization. 



keT 


\Ef\ 

Figure 4. The Temperature dependence of xvt for positive Fermi energy (red curve) 
and negative Fermi energy (blue curve). 


approaches a constant value 1. In the intermediate regime 
monotonically decreasing as Ep/ksT increases. 

Eq. fl37|l can be re-arranged as 

Tr‘^2\EF\ 


-5 < Ep/ksT < 5, Xvt is 


a, 


Xvt — To—2NokB- 
n 3 


E 


R 


1 - 


b{t2) -t2a{t2) 


(38) 


To make the temperature dependence of xvt clear, we replot xvt (in units of 
TopNokB^^^) as function of ^ (^) in Fig. 4. 

It shows that xvt increases with increasing temperature. When kBT/Ep < 0.2, 
Xvt fully vanishes, when kBT/Ep continues to increase the contribution of the band 
valley structure will be included and dominates Xvt, then xvt is enhanced. When 
kpT/ {—Ep) < 0.2, XvT/vof is almost exactly linear in kpT/ {—Ep) 

with slope 1. When the temperature increases, not only the band valley contributes, 
but also the electron states above the band crossing point are included due to the 
thermal smearing, thus the TISP is suppressed. Moreover, when 3> 1, ^2 —t 0, 


) b{t2)-t2a{t2) 


-)■ 


in T. 


so Xvt > 1 


\Ef\ 

^ Xvt{Ep = 0) = TopNokB^h^ linear 


3 2 E 


R 
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In conclusion, we have calculated the thermoelectric responses of spin polarization in 
2D Rashba system. By self-consistently determining the transport time, we exactly 
solved the Boltzmann equation when static impurity scatterings dominate the electron 
relaxation process. It was shown that the electric-held induced spin polarization is 
linearly dependent on the Fermi energy when only the lower band is occupied, different 
from the Edelstein behavior when both bands occupied. Higher spin polarization 
efficiency is achieved when the Fermi energy lies below the band crossing point. It was 
found that the temperature-gradient induced spin polarization continuously increase to 
a saturation value as the Fermi energy decreases below the band crossing point. In 
addition, the temperature-gradient induced spin polarization tends to zero at vanishing 
temperatures. 

This work may stimulate more experimental and theoretical works on the electrical 
and thermal spin control in Rashba semiconductors BiTeX (X=C1, Br, I) and BiTeX 
quantum wells, as well as other materials with giant Rashba spin splitting. 
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